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Abstract: We examine the moduli dynamics of a specific class of super gravity-inspired 
BPS braneworlds, clarifying the role of bulk scalar fields in brane collisions. The model 
contains as a special case the Randall-Sundrum model both with and without a free, 
massless bulk scalar field. Its low-energy effective theory is derived with a moduli space 
approximation (MSA) and agrees with the corresponding results derived elsewhere. Rather 
than stabilising the radion, we look at cosmological evolution of the system stimulated by 
breaking the BPS condition on the branes. We examine in detail the range of validity of 
the MSA in both the RS and BPS case, paying particular attention to the divergences that 
can arise during a collision of the branes. In the absence of perturbations such an event 
is finite in the RS model, and accurately described by the low-energy effective theory. We 
demonstrate, however, that a collision is divergent in the BPS case even with an exact 
FRW geometry. 
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1. Introduction 

Braneworld models for the Universe have recently been the subject of considerable theoret- 
ical interest [1]. Much of this research has been inspired by the work of Horava and Witten 
[2], who showed that the strong coupling limit of E8 x E8 heterotic String Theory could 
be described by an eleven-dimensional supergravity with the eleventh dimension made up 
of the orbifold S 1 /Z2, i.e. an interval with reflection symmetry about its endpoints. These 
endpoints define ten-dimensional submanifolds, the branes, which sit at the orbifold fixed 
points and define the boundaries of the spacetime. The other six dimensions could consis- 
tently be compactified on a Calabi-Yau threefold whose characteristic scale is considerably 
smaller than the interbrane distance. This theory then led to many toy models of spacetime 
as a five-dimensional manifold (the bulk) bounded by two branes with a Z2 symmetry. 
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One of the most well-known of these is the Randall-Sundrum I model [3]. Whereas the 
Horava-Witten model allows a large number of fields to propagate in the bulk, this model 
has only a bulk cosmological constant. The two boundary branes carry equal and opposite 
tensions whose magnitudes can be fine tuned such that the effective cosmological constants 
on the branes vanish, provided the bulk cosmological constant is negative. The resulting 
low-energy theory can be derived in a number of ways [4, 5, 6] and can be formulated in 
terms of a single modulus field, the radion, representing the proper distance between the 
branes. This four-dimensional low-energy theory describes the physics as viewed by an 
observer confined to one of the branes and, in the case of exact cosmological symmetry, 
turns out to be perfectly well defined even when the size of the fifth dimension vanishes, 
i.e. the branes collide. 

The radion field is classically massless, representing the fact that the brane positions are 
arbitrary when the brane tensions are fine tuned. This is phenomenologically undesirable, 
since a massless scalar field is not observed in nature. By detuning the brane tensions 
from their preferred values one can generate a potential for the radion; one could imagine 
these detuning potentials on the branes coming either from a more consistent quantum 
mechanical treatment of the scalar field [7] or from sort of SUSY-breaking mechanism. In 
this paper we shall just insert it by hand. 

We shall consider a supergravity-motivated generalisation [8, 4] of the Randall-Sundrum 
model where a scalar field \P is allowed to propagate in the bulk. Its potential U replaces 
the Randall-Sundrum bulk cosmological constant and the 'superpotential' induced on the 
branes V replaces the tensions. Again the branes can be Minkowski if V and U are re- 
lated by a BPS condition, the analogue of the Randall-Sundrum fine tuning. In fact, this 
model reduces exactly to the Randall-Sundrum model if its free parameter a, defined by 
V oc exp a^>, vanishes. This time, the four-dimensional effective theory contains two moduli 
fields and, by detuning the brane potentials, one can generate dynamics and cosmological 
solutions as before. 

The cosmological consequences of this model have been explored in detail in [4]. In this 
work we focus on brane collisions, events which have been used recently in the literature 
to try to illuminate and resolve the Big Bang singularity [9, 10]. It is already known that 
perturbations diverge during such an event; in this paper, we investigate how the simple 
generalisation to BPS bulk scalar fields effects both the dynamics of the radion and the 
regularity of the collision. Such bulk scalar fields arise naturally in supersymmetric theories 
[8]. 

This paper is organised as follows. In §|| we describe the model and derive the exact 
projected Einstein equations on the brane. In the cosmological context we obtain the 
Friedmann equation for the induced FRW brane geometry. This is local (i.e. defined purely 
in terms of quantities on the brane) apart from the well-known Weyl tensor representing the 
gravitational effects of the bulk. Discussing the low-energy effective theory in §|3] connects 
this quantity to the moduli fields and allows one to write down an approximate, closed set 
of equations for the motion of the branes. The dynamics in the Randall-Sundrum case are 
reviewed in §||, with the identification of the Weyl tensor explaining the finiteness of the 
system through a brane collision. The more general case a ^ is studied in §|5|, showing 
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that in all but the most contrived case the scalar field will diverge during a brane collision. 
The conclusions are summarised in §||. 

2. The Model 

A general action for the two-brane system is 



S 
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ignoring for the time being the presence of possible matter actions for the branes. and 



h 



(2) 



ab are the induced metrics on the positive- and negative-tension branes respectively, U 



and V > are the bulk and brane potentials and K is the trace of the extrinsic curvature 
K ab taking outward-pointing normals. The factor of k 2 ensures that $ is dimensionless. 

The index convention to be used is that the metric signature is 'mostly plus' — h + ++, 
a,b,c... are five-dimensional indices and fi,u,p... are four-dimensional, running from to 
3. The unusual factor of 2 in the Gibbons-Hawking sectors is because of the Z2 orbifold 
nature of the bulk, and the five-dimensional integral in the action is taken to mean an 
integral over two copies of the bulk spacetime Region I between the two branes, see Fig.|l[ 
As described in the introduction the bulk geometry is assumed to be Z2-symmetrical about 
the position of the two branes, which simplifies the junction conditions. The 1 and 2 on 
the brane integrals refer to positive and negative tension respectively. 

Variation of the action with respect to the metric and the field yield the usual bulk 
equations and junction conditions for each brane [11]. Considering the positive-tension 
brane (the analysis for the other brane is equivalent with V - 
and find 



-V), we define h ab = h$ 



G ab = \d a <S>d b <5> - - A g ab d<$> 2 - \u{<b)g ab 
<i<3? 



(2.2) 



[K ab ] = -4 T ab - -Th ab 



[n-m = 24^- 



where 



T ab = -V($)h, 



ab 



(2-3) 



and [X] means the jump in X in the same direction as the normal (hence both these 
quantities are invariant under n — > —n). Assuming Z2 symmetry the junction conditions 
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become 

- K ab = -l4v(<S>)h ab (2.4) 

- n -d$ = 4^- (2.5) 

where quantities are evaluated at the edge of Region I, which the normal points outward 
from as depicted in Fig. [I]. The 4D Einstein tensor G ab (h) is given by the standard result 

G ab = \ \G cd h c a h d b + (G cd n c n d - ]g] h ab ) (2.6) 



+KK ab - K c a K bc ~\{k 2 - K cd K cd ) h ab - E ab 

where n a is the spacelike unit normal to the brane and E ab is the (traceless) electric part 
of the Weyl tensor 

E a b = C a cbdn C n d 




Figure 1: Bulk and boundary structure: The bulk topology is E x S1/Z2, with the branes 
sitting at the orbifold fixed points. There are therefore two identical copies of the bulk, Regions I 
& II, with four boundary planes in all as shown, giving rise to two copies of the bulk action and of 
the GH boundary terms. There is only one copy of each brane worldvolume (and matter) action. 

From Q, 

G cd h a h d + (G cd n c n d - hs\ h ab 

= Jv a $V fe $ - ^V^ 2 h ab ~\~ —— (71 d$) 2 h ab - \uh ab 
2 lb lb 8 

where V is the covariant derivative of the induced metric h, and from (|2.4j ), 
KK ab - K c a K bc -\{K 2 - K cd K cd ) h ab = -^4v 2 h ab 
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Substituting all this back into fl2.6|) with the junction condition ( j2.5| ) we obtain the pro- 
jected Einstein equation on the brane: 




(2.7) 



If matter were included on the brane via 



V($)h ab + r ab 



(2.8) 



then there would also be the terms 



(2.9) 



where 7r a f, depends quadratically on T ab . This has two consequences; in order to recover the 
Friedmann equation we must live on the positive-tension brane (at least at the classical 
level), and be at an energy scale much less than the brane tension so that the resulting 
Tq term can be regarded as a small correction important only in the early Universe. Note 
that the effective four-dimensional cosmological 'constant' is given by 



Such self-tuned potentials often arise in the context of supergravity models where U and 

V would be derived from the same superpotential [8]. We shall for convenience refer to 
such potentials as being 'supersymmetric' and the function V as the superpotential. This 
relation between the bulk and brane potentials is a generalisation of the usual Randall- 
Sundrum (henceforth RS) fine tuning, which can be reproduced by taking U = — 2A and 

V = a; ( |2.10| ) then reduces to the familiar form 



3. Low-energy Effective Theory 

Whilst ( [2.7D gives the exact Einstein equations for an observer on the brane it has little 
predictive power, since it contains the term E ab which is not defined in terms of data on the 
brane. Such explicit dependence on the bulk geometry is also present in the Klein-Gordon 
equation ( |2.3D which, when written out explicitly in terms of four-dimensional covariant 




which will vanish for potentials satisfying 




(2.10) 
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derivatives, will contain a term of the form (n • d) 2 <£, not prescribed by the value of the 
field on, and its derivatives along, the brane. 

For this reason one seeks a four-dimensional effective theory in which these non-local 
quantities are replaced by scalar 'moduli fields'. Such low-energy effective theories can 
be derived in many different, but essentially equivalent ways. Here we shall employ the 
moduli space approximation; the result ( 3.14 ) agrees with that obtained by perturbative 
expansion [12]. 

Note that a four-dimensional description, low energy or otherwise, in terms of an 
effective action cannot hope to reproduce a Friedmann equation with the terms quadratic 
in the brane energy-momentum tensor given by ( |2.9| ). For example [13, 14], in the RS case, 
these can be viewed as arising approximately from the trace anomaly in the CFT defined 
on the branes in the context of the AdS-CFT correspondence, and hence cannot be derived 
from the variation of an action. Higher order derivative terms could be included in the 
effective action; these will approximate then more and more closely the effects of E^ v in 
the projected Einstein equations, but the quadratic stress-energy terms cannot be obtained 
this way. 



3.1 BPS backgrounds 

The significance of supersymmetric (henceforth SUSY) potentials is that they allow privi- 
leged configurations of the system where the brane positions are arbitrary. These configu- 
rations are, in a sense, a ground state of the system, with a high degree of symmetry. We 
therefore look for solutions for the metric and the scalar field which do not depend on the 
transverse directions (i.e. are static and transversely homogeneous). We take Gaussian 
Normal coordinates away from the positive-tension brane 

ds 2 = a{z) 2 r]^ u dx^dx v + dz 2 

where z is the proper distance along the normal to the brane (increasing towards the other 
brane) and x 11 parameterise the flat, transverse foliations. Since we are assuming fine- 
tuning of the brane tensions the transverse foliations are flat, hence the use of r\^ v above. 
Due to the symmetry of this 'vacuum' configuration the branes are at constant z. In the 
next section, when we relax the requirements of staticity and homogeneity, we will consider 
more general coordinate systems. Taking <$(i, x, z) = <&(z), the Einstein and Klein-Gordon 
equations become: 

-r + — = -t^$ - t:U 



a 2 a 12 6 
a 2 24 12 



= XJ (3.1) 



+ 4— $' = — 
a d<3? 
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where ' = d/dz. For potentials U satisfying ( |2.10 ) this can be written as 

1 



a 

a r 
i\ 2 

a 
a 



6 



12 



24 24 " 



1 

36 



(3.2) 
(3.3) 



noting that there are only two independent equations in (3.1). An interesting family of 
solutions follow from the first-order system 



a 
a 



4v 



,dV 



(3.4) 



6 ' ~ °d$ 

which can easily be seen to solve the full Einstein-Klein-Gordon equations fl3.2j]3.3j ). In the 
underlying supergravity theory [8] these are BPS configurations; the 'no force' condition 
between BPS branes manifests itself in that ( |3.4 ) implies the junction conditions, which 
are given by 



a 
a 



4v 



6 



*/| 2 dV 



(3.5) 



where the quantities are evaluated at the positive-tension brane and we have used K a b = 
g' ab /2 for surfaces of constant z. In other words, there exist static solutions of the system 
for arbitrary brane positions. A solution of ( \5-4 ) will from now on be referred to as a 
'BPS background'. Note that the free parameters in these solutions are the value of the 
scalefactor a and the scalar field on the positive-tension brane. 

Note that, in the metric junction condition ( |3.5P , a' /a is proportional to the extrin- 
sic curvature on the brane. Since the normal to the negative-tension brane must point 
in the opposite direction for consistency (i.e. either both outward or both inward), the 
corresponding junction condition there will pick up a minus sign. Therefore, in order for 
the junction conditions to be satisfied at both branes simultaneously, the negative-tension 
brane must have tension V-(&) = —V+($>) as is assumed from the start in the action (|2.1|), 

As we shall see in the next section, it is desirable to detune the brane potentials from 
their 'supersymmetric' values. This will generate a potential for the moduli fields in the 
four-dimensional effective theory, which is desirable for both non-trivial dynamics and some 
hope of their stabilisation. We shall insert this by hand, 



Vz($) = -V(<f>)+w($) 



(3.6) 



where V\ 2 are the tensions on the two branes, V is the supersymmetric value given by 



(|2.10[) and v, w are small perturbations, v, w -C V. For definiteness, we shall consider only 
exponential superpotentials, 

(3.7) 



V = % a<s> 



which reduce to the RS model on taking the limit a — > whereupon the constant k, which 
has dimensions of inverse length, is the curvature scale of the then AdS bulk. 
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3.2 The Moduli Space Approximation 



In the previous subsection we derived a set of equations (3.4) whose solutions we will use as 



a ground state for our model. The corresponding metric and scalar field profiles were static 
and homogeneous, depending only on the bulk coordinate. This profile for the scalar field 
and the metric satisfies the junction conditions ( |3.5[) at the branes provided the branes are 
parallel and static. In order to examine small perturbations around this vacuum we allow 
the brane positions to fluctuate as depicted in Fig. [T|, and incorporate the graviton zero 
mode [8] by generalising the metric ansatz to 



ds 2 = dz 2 + a(z) 2 g^(x)dx^dx u (3.8) 



In this coordinate system the perturbed brane positions are given by z = z\{x) and z = 
Z2{x) (from now on, subscripts 1 and 2 will refer to evaluation at the positive- and negative- 
tension branes respectively). The radion, the proper distance between the two branes, is 
given by r(x) = z<i — z\. This is the local size of the fifth-dimension; the only physical 
spacetime is that between the two branes. 

We shall assume that any ^-dependence is small, i.e. that transverse derivatives are 
much smaller than normal derivatives. Also, we shall assume that the SUSY-breaking 
potentials v and w (and matter actions were we to be considering them) are small, to be 
consistent with the fact that we must necessarily lose the quadratic terms in the Fried- 
mann equation in using an effective action, as discussed above. In this sense then the 
effective four-dimensional theory we will obtain is only valid at low energies, specifically 
for VjWjTfj,!, <C V. Radion fluctuations can be regarded as a small perturbation around the 
BPS configuration; by substituting this perturbed brane positions back into the action, 
keeping the BPS profile for the metric and scalar field, we shall obtain an effective theory 
governing the motion of the moduli fields Z\ and Z2 ■ 



For the superpotential (|3.7|) the BPS profile is given by 

a(z) = C(z) 1/6a2 

$(z) = --log* 
a 

£{z) = [6ka 2 (zq - z)] 



= --log£(*) (3.9) 
a 



where zq is a constant of integration and we have assumed As is generally the case 

with self-tuned potentials there is a singularity in the bulk at z = Zq. The theory will break 
down when this singularity lies in the physical region of spacetime between the two branes 
[4]- 
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The action is given by 



Sfuii = S g + S$, + Sb + 5gh 

5g = i / d " x ^R{g) 



-d$> 2 - 



S B = J d*xV-hW (-V - v) 
+ J^d 4 xV-h( 2 ) (v-u?J 
Sqh = —7) [ <Pxy/-hP-)K 1 + -i f d 4 xV-hWK 2 



(3.10) 
(3.11) 
(3.12) 



(3.13) 



remembering that 



d 4 x 



z\(x) 



dz^J —g a 4 



Keeping the BPS profile for the metric and the field in place, we can evaluate this term 
by term to obtain a four-dimensional effective theory in terms of the Z{. For v = w = 
all terms not involving x-derivatives will sum to zero since they would otherwise generate 
a potential term for the moduli fields or a cosmological constant term, which both must 
vanish in this fine-tuned case. Hence we can discard any terms not involving transverse 
derivatives. 

The result, expressed in terms of the induced metric on the positive-tension brane 
h^u = hn U , is derived in the Appendix, and is given by 

S = 777—2 / d 4 xV^h [n 2 R(h) - 7 AB VU • V£b - V] 



with 



& = 



1 + 3q 2 



l/3a 2 



7 



,11 



7 



12 



[6a 4 (l + 3a 2 )] 
— [6a 



-i 1 



41-1 1 (& 



7 



22 



= M 



-i 1 



I) 



a l 

l/3a 2 



3a 2 + 



l+l/3a 2 



(3.14) 



l/3a 2 



V = 2k4 



v + w ( — 
?i 



C \ 2/3a 2 ' 
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in agreement with the result of [12] (using different normalisations for the fields). Here, and 
elsewhere, covariant derivatives are taken with respect to the metric h^ u unless otherwise 
specified. Finally we put the action into a more easily interpretable form, in terms of the 
radion and value of the scalar field on the positive-tension brane. For convenience we make 
the connection to the radion via the approximate conformal factor uj 2 relating the two 
induced metrics, 

w2= /a^y (JU5) 



and define 



ai(x) 



#r) = l-^ 1 + 3a2 ) = l- f|j (3.16) 

so that the branes coincide for ip = whilst ip = 1 signifies an infinite redshift between the 
two branes; this could either mean their proper separation is infinite or that the second 
brane has hit the singularity. The value of the scalar field on the positive-tension brane is 

»7(x) = --log£i (3.17) 

a 



l+l/3a 2 



giving our final action as 

.<? = 

2k 2 



S = At / d l .rv'-hc-' ul 



ipR(h) 



2(3 (1 - V) 

- J d 4 xV^h [v+(l- ^) 2//3 u>} (3.18) 

where 

K 2 = kK\(3 (3 = l + 3a 2 (3.19) 

Since the underlying action has five-dimensional coordinate invariance, the choice of co- 
ordinate system used to perform the above dimensional reduction cannot affect the final 
action, which is a scalar quantity. Furthermore, if one changes the coordinate system so 
that neither, one, or both of the branes are at fixed positions, each component of the action 
is still separately invariant for the same reason. This will, of course, depend on which four- 
dimensional metric is being used, but if one is careful to express everything in terms of the 
same metric (h^ u , for example) then it is easy to check that each contribution to the total 
effective action is unchanged. For example, one could use a coordinate system in which 
the branes are fixed, say at y = and y = 1, via the transformation z = z±(x) + r(x)y. 
The scalar field modulus then enters as usual, being defined as the value of the field at, for 
example, y = 0, but this time the radion enters directly in the metric, 



ds 2 = r{xfdy 2 + 2yd ll r{x)dx ll dy (3.20) 



+ 



a{z(x,y)) 
aAx) 



2 



h^ u (x) +y df,r(x)d u r(x) 



dx^dx u 
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The presence of these extra terms in the metric ensures that the boundary terms still give 
precisely the same contribution. 



3.3 The RS Limit 



From (3.9), we see that 



a(z) ~ aoe 



7] = $(Zl(x)) — ► 



(3.21) 



as a — > 0, i.e. one recovers a AdS profile in the bulk and the scalar field vanishes. Hence 
the RS model should be understood as being recoverable at the level of the action by taking 
a = r\ = 0: 



which is the standard action found in the literature [5]. Note that setting ij = is 
consistent with the equation of motion obtained from the full action in the limit a — > 0: 



If one takes a = does not set r\ to zero, the system describes a RS braneworld with 
a free, massless bulk scalar field added by hand. With a different normalisation for 77, the 
resulting action is in agreement with the result of [6]. 

3.4 Scalar Degrees of Freedom as Goldstone Bosons 

The moduli fields have an interesting interpretation in terms of symmetry breaking. In the 
RS case, the radion can be thought of as the Goldstone boson associated with the breaking 
of translation invariance. For example, if we fix the coordinate gauge invariance (in the bulk 
direction) by defining the positive-tension brane to be at z = 0, there is still a continuous 
family of 'vacuum states' parameterised by the arbitrary position of the negative-tension 
brane at z = z\. Choosing a particular value of z\ breaks this symmetry, giving rise to 
a Goldstone mode which can be identified with the radion. In the presence of the scalar 
field, one is also free to choose the value of the scalar field e.g. at z = 0, $xj there are then 
two continuously deformable parameters z\ and <l?i parameterising the 'vacuum manifold', 
giving rise to two massless degrees of freedom when a specific choice of ground state is 
made. When the tensions are detuned the moduli develop a potential, and hence a mass, 
since they are no longer Goldstone bosons - there is no longer a continuous family of static 
solutions due to the violation of the junction conditions. 

4. Dynamics for RS 

In this section we re-derive the well-known result that tension perturbations in RS generate 
a potential for the radion (see, for example, [10, 19]) since it will be useful for the next 
Section where we allow 




(3.22) 



i\)Ut) + Vtp ■ V?? = 



(3.23) 
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4.1 Equations of motion 

The equations of motion which follow from the RS effective action for the positive-tension 
brane (which we shall be working with for now on unless otherwise specified) are 



iPG 



fJtV 



+ 



3 



2(i-V0 



k 2 h 



V + (1 — ifi) W 



Dip 



1 Vip 2 



4^ 

2i-v ~ ~ 



(1-V) (v + (l-*l>)iv) 



(4.1) 



(4.2) 



where now ip is given by tp = 1 — exp(— 2kr). K,as defined from the action, is given by 

(4.3) 



which is also the effective gravitational coupling on the brane; this can be identified directly 
from the projected Einstein equations in the presence of matter (2.7,2.9), which give 



^4 



-V 



(4.4) 



For the RS case, V = 6k /k 2 , this is consistent with (|4.3|) , i.e. k = k±. In the absence of 
matter, the projected Einstein equations for the positive-tension brane are 



G 



(4.5) 



allow the identification, to this level of approximation, of 



E. 



fJ,U 



h 



Dip 



1, 



v + (1 — ip) w 



2, 

-t- K n,fj, u —T— 



(4.6) 



from ([4.1]). Note that the tracelessness of E^ u implies the equation of motion ( |4.2j ) for ^. 
In fact [13], in the absence of matter, the requirement that a single scalar field coupled 
to gravity yields Einstein equations consistent with the tracelessness of the Weyl tensor 
determines that, defining the scalar field to be the coefficient of R as above, the action 
must be of the form 



S oc / d x 



-h 



ipR 



2(1-V) 



Vip 1 



+ A + B(i-ipy 



(4.7) 
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When the system is fine-tuned (A = B = 0) one can write the resulting Einstein equations 
as 

G„ 



2rrflj) 



and comparing with ( |4.5| ) one can identify the effective energy-momentum tensor for i/j 
with —Euv The fact then that T^ v must be traceless is sign of the underlying conformal 
invariance of the RS effective action [14]. 

4.2 Dynamics 



® 




X 



Figure 2: Effective potential for radion, RS case: The dashed plot has v = 0.1, w = 0.2, and 
the potential drives the branes to infinite separation \ — 0. The solid plot has v — 0.1, w — —0.2, 
giving an unstable fixed point at % c = 1/ y/2. The branes either collide or move apart depending 
on the initial conditions, k — 1 in both cases. 

The field ip has a non-standard kinetic term and, hence, its dynamics are not immedi- 
ately obvious from ( [4 . 2[) . Defining 



X 



- kr 



we find 



_ 4k 2 

□X = —X {v + X 2 w) 



i.e. the field \ moves in the potential 

Hx) = 

This has a turning point at 

v i 



4k 2 /l 2 1 4 
— {2 XV+ 4 XW 



xl 



w 



l+ V - 

w 



(4.8) 
(4.9) 

(4.10) 

(4.11) 



which is only inside the physical range < \ < 1 for v and w of opposite sign and \ v\ < \w\. 
This stationary point is unstable for a de-Sitter positive-tension brane (i.e. positive effective 
cosmological constant, v > 0), i.e. an inflating brane, and stable only when the brane is 
collapsing (v < 0). Depending on the initial conditions the branes will either be driven 
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apart tox = 044>r^ooor will collide, x = 1 ^ r = 0. However, if there is no stationary 
point in the potential (for example, if v > 0, when v + w > 0) the branes will always be 
driven apart. The potential ( 4.10 ) is plotted in FigJ^ for two different cases, and FigsJ^ 



and [| show explicitly the above behaviour. 
4.3 Brane Collisions 

As we have seen it is possible to produce an initial configuration which will lead to a brane 
collision, which corresponds to ip = 0. Whilst this appears in general to be a singular 
point of the equations of motion, in the special case of an exact FRW induced geometry 
the collision (in the Brane Frame, as used here) is regular. We take a FRW metric with 
scalefactor a(t), curvature K = 0, ±1 and Hubble parameter H = a/a. The Bianchi 
Identity, 

VG^ = => VE^ = 0^E 00 k a" 4 (4.12) 



H 2 + K/a 2 = — + —v (4.13) 



gives, from (|4.5|) , 

C_ K 2 

for some constant dark radiation coefficient C. Hence the Hubble parameter does not 
diverge in the collision; the equations of motion ensure that, despite appearing to be 
singular at ip = from ( |4.6p , Eqq actually only behaves as dark radiation. This is why we 
have chosen to work in the Brane Frame, rather than the Einstein Frame defined by 



hfiu = tphfiu ds 2 = tp [dt 2 — a 2 dx 2 ) 

which gives 



= dP — d 2 dx. 2 , 



a dt v¥ 2^ 3 / 2 
which will, in general, diverge as rp — ► 0. 

It is therefore possible to follow the evolution of the system through a brane collision 
by using (or, rather, its derivative) in place of ( |4.1| ). Taking K = for simplicity, the 
equations of motion for ip and H are 

H = -2H 2 + — v (4.15) 
3 

$ = -SHip ' ' " 



21-ip 

4k 2 / \ 

+ _(1-^) (v + (l-ip)w) (4.16) 

which are manifestly regular as ip — > 0. Note that, for v > 0, ( 4.15| ) has analytic solution 



H = k J | tanh ( 2k J |f + cnst ) ( LIT) 
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i.e. the Universe rapidly approaches de Sitter space, H —> = Ky/v/3. The evolution of 
the Hubble parameter is completely independent of that of the radion and, in particular, 
is not affected by tp = 0. 

The initial value of H can be determined from the 00-component of (4.1), which gives 



H 



+ 



■ (W(i-v) 2 ^) 



(4.18) 



which in general has two solutions. Some possible evolutions of H and ip for different initial 
conditions are given in Figs.^ and ||] which uses the potential of Fig.|| (ip c = 0.5, H*/k 
0.183). (|4l5|) can also be obtained as a linear combination of the three equations of motion 
obtained from (O) and Q4.2| ); using these equations, ( |4.16| ), ( 4.18 ) and 



2H 2 + 3H 2 = H 



± 

fjt 4^(1 - V) 

2 



(4.19) 



gives the same results but more care must be taken with the numerical integration to 
cope with the divergence of some terms as if) — > 0. 




Figure 3: Brane collision in RS: ip(t) and H(t) are plotted for initial conditions leading to a 
brane collision (-0 = 0). H(t) — > iJ* as expected. 



4.4 Accuracy of the MSA 

As we have seen, the four-dimensional effective theory contains a single modulus, ip, related 
to the local radius of the orbifold dimension. The effective theory assumed that Vzf/a 2 <C 
1, and we need to check that these conditions are not violated during the evolution of the 
system. It is actually sufficient to check that Vr 2 /w 2 <C 1- Taking a Gaussian Normal 
coordinate system where the positive-tension brane lies at y = 0, the metric is given by 

ds 2 = dy 2 + b(y) 2 h^dx»dx v (4.20) 

where the scalefactor b{y) is normalised to 6(0) = 1, hence hu U is indeed the induced metric. 
The same effective action ( |3.22| ) is then obtained provided Vr 2 /6 2 <C 1, where b(x) is the 
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t 



2 4 6 8 10 12 14 

Figure 4: Branes driven apart: ip(t) — > 1 asymptotically for suitable choice of initial conditions. 
Note that e ceases to be small as ^ — » 1 



scalefactor at the second brane, y = r(x). Clearly we can identify = u(x), giving the 
above condition on Vr 2 . In other words, although only checking Vr 2 to be suitably small 
ignores the possibility that z\ and Z2 could oscillate wildly but coherently (i.e. that Vr 2 
small =£• Vz 2 small), all that is of physical relevance is the separation, r(x). In terms of 
the modulus tp, the condition is 



|v-W| = w 



Vt/> 2 



< i 



(4.21) 



As can be seen from Fig. |], this suggests that the MSA will become inaccurate as the 
branes move apart, t/j — > 1; a given fluctuation will cease to be small as lengthscales shrink 
on the negative-tension brane. As tp — > though, we just need that \7ip 2 <C k 2 ; the 
approximation appears to be valid through the collision. In all the numerical simulations 
we shall define our units by k = k 2 = 1. 



5. Dynamics with Bulk Scalar 

Since the scalar field is assumed to take its BPS profile in the background the effect of 
non-zero a on the dynamics is just to alter some of the coefficients in the potential. In 
particular, we are not attempting to produce a potential capable of stabilising the radion 
as in [18, 19]. The extra degree of freedom, however, gives rise to key differences in the 
brane collision. 

The gravitational coupling of matter on the brane is now given by 

k 2 M = kK 2 e^ = je a * (5.1) 

For simplicity we shall consider still the case where the tension perturbations v and w 
are constants; qualitatively identical results are obtained for perturbations with the same 
functional form as the superpotential, i.e. v(<&) = 5vexpa<&. 
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5.1 Equations of Motion 



The variation of ( 3.1S| ) with respect to h^, ip and rj gives 



R 



tpR 



VfjVu^ - h^Uip - aip (V^VuV - h^Urj) 
-2a (V^VV^ry - h^Vip • Vr?) 



+ 



1 



1 



2/?(l-# + ^ + 



2/3 (1 - VO 
2 

rcV^V (?; + (l-^) 2//3 u;) 

3 V^ 2 1 2 3aV??-VV> 
' — Vrr 4- 
2 

T 



3 



/3 (1 - V) 
3 VV> 2 



(5 (1-V) 
e a?? (1 - V) 2//3_1 



2 a a 



(5.2) 



respectively (/3 = 1 + 3a 2 ). Eliminating i? from the last two of these equations yields 



dip = aS/r] ■ Vip 



vv> 2 



2(1-V) 



Ak 1 



(l-ip)v + (l-ip) 2//3 w 



Or/ = aVry 

4a 
1 



2 Vt?-VV> 3a V^ 2 



V> 2/3 



+— « 2 e a?? u 



(5.3) 



(5.4) 



consistent with ( [4.1| ) and ( |4.2| ). The behaviour of ^ is qualitatively similar to the RS case, 
with the critical value of tp modified to 



l-3a 2 



(5.5) 



which is independent of r/. Although there is no corresponding stationary value of rj unless 
a = 0, ( |5.3|) ensures ip = ip c , Vtp = =^ Dip = 0. Note that (|5.3| ) is regular as tp —> 0. 

The influence of the bulk scalar field is therefore not great on the effective potential 
in which the radion moves; however, we shall see that it has a more important role to play 
during a collision. 
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5.2 Brane Collisions and Cosmological Evolution 

Taking again a flat FRW metric, the Einstein equations ( |5.2| ) give the Hubble constraint 



H 1 



and the evolution equation 



2H 



ih 1 th 2 

-H^ + aHr] + — —^ 
ip 4/3 tb(l-ib) 



ft 2 
+ I2 



(5.6) 



3ib 



-3F 2 - t 



a 2 + -)v 2 



(v + (l-ibf /l3 w 



tb 

2H^- + aij + 2aHf] 



+ 2<M 
4/3ib(l-ib) ib 



-« 



3H Z + H 

2 



(y + (1 - V) 



2//3 



3-0 



3 m? 



aHf] 



w 

r) 2 
T 



1 



V 2 



4/3-0(1 -V0 



| )« + (! 



using (5.3). Finally, substituting in from (5^6) gives the analogue of ( 4.15 ): 



H = -2H 2 -£ + %e^v 
12 3/3 



(5.7) 



(5.8) 



So again we have manifestly finite equations of motion for H and tp as ij) — * 0. However, 

— > oo then, from (|5.S|), there 



the equation of motion (|5.4| ) is not free of divergences; if |r/| 
will be a corresponding divergence in H. Such behaviour is shown in Fig|| 
For the a = case, i.e. RS with a free, massless scalar field, (5^) gives 

1 



V M (V'V^??) = 



77 oc 



(5.9) 



Hence 77 will tend to either a constant (if 77 = 0, since the field 77 then has no time dependence 
and decouples) or to ±00 as ip — > depending on the initial sign of f\. 

In the general case it is difficult to make much progress analytically. The numerical 
results for four sets of initial conditions and values of v, w and a are shown in Fig|| 
appearing to show that 77 diverges at the collision with a sign opposite to that of a (a 
conclusion supported by thorough numerical investigation, of which Fig.|| is only a sample). 
This implies that the tension on the branes, given by ( |3.7| ), vanishes as the branes collide. 

The conclusion is then that, although the scalar field does not have much impact on 
the dynamics of the radion itself, it will, in general, diverge during a collision (one could 
envisage a situation where it would not, for example, by setting a = f/ = 0, in which case 
77 would just remain constant and the collision would be regular). This is similar to the 
evolution of perturbations, which are found to diverge logarithmically in the RS case [10], 
a feature common to scalar-tensor theories when the gravitational constant changes sign. 
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When the initial conditions are such that the branes do not collide, ip is rapidly driven 
to 1. 7] is then approximately governed (for a 7^ 0) by 



rj + 3Hr] ~ afj 2 —K 2 e art v 



ft 



(5.10) 



The Hubble constant is slowly driven to zero whilst \rj\ grows approximately logarithmically 
(a full numerical solution is given in Fig]| for a = -0.2). There is no analytic solution 
to the system even if one approximates ip = 1; one can, however, predict the late time 
behaviour 



a 



■log|l + t-to| 



(5.11) 



with the timescale for the transition to the logarithmic behaviour increasing without limit 
as a —* 0. This is clearly demonstrated in Fig. ^. 

5.3 A Caveat 

As before, the self-consistency of the effective theory can be checked by monitoring the size 
of the quantities 



f'2 



Vz 2 



,(l/3a-2a)v 

36k 2 a 2 
,(l/3a-2a)r? 



|Vr? 2 



36k 2 a 2 



(1-V) 



-3/(l+3a 2 ) 



x| (1 - V) Vrj + 3aVij\ 2 
~ ei as ip — > 

However, due to the divergence of 77, these are not necessarily small during the collision, 
as can be seen in Fig.||. The reason for this divergence, and the reason for its absence in 
the RS case, can be understood as follows. In our coordinate system in which both the 
bulk and the scalar field are static, one can interpret a divergence in 77 as a divergence, or 
a tending to the singularity depending on the sign, of the coordinate z of the positive- 
tension brane. In the RS case, the two coordinates diverge in the collision. However, the 
effective theory can be expressed solely in terms of the difference between the two values, 
the radion, which remains finite. In the more general case, the divergence of the second 
modulus cannot be removed in such a way. The breakdown of the effective theory can also 
be understood by rederiving the effective action in a Gaussian-Normal coordinate system 
about the positive-tension brane, analogously to the procedure used in [6]. One performs a 
perturbative expansion about a static ground state [12], requiring transverse derivatives of 
the scalar field to be small. However, this procedure will break down when this condition 
is no longer met, in particular when 77 diverges during the collision. 

Whilst then we cannot trust the MSA right at the moment of collision, we have a firm 
numerical tool for identifying dynamically the region in which it is, indeed, valid. We can 
use the MSA to identify the approach to a collision and, although we have not ruled out the 
possibility of higher-order effects repelling the branes again, it seems likely that a collision 
will then take place, accompanied by a divergence in the bulk scalar field. 
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5.4 Einstein vs. Brane Frame 



In this paper we have chosen to work in the Brane Frame, defined by the induced metric 
on the positive-tension brane, since this has the best regularity properties at the brane 
collision. For late-time evolution the Einstein Frame is most often used in the literature; 
from the action ( |3 . 1 8| ) , this can be read off as 

= e~ av W ( 5 - 12 ) 

As V —> the two metrics will give wildly different physics, but in the RS case, with or 
without the scalar field, the two metrics become identical as V — * lj s ° the distinction 
between the two frames becomes unimportant. However, for a / 0, we see from (5.11) 
that 

e -m V oc 4 as ip -»■ 1 (5.13) 

so that, even for late times, there is still a distinction between the two frames. The 
cosmological evolution of the moduli fields in the Einstein Frame are investigated in detail 
in [4]. 



6. Conclusions 



In this paper, we studied the evolution of the moduli fields of a particular class of SUGRA- 
inspired braneworld models, containing the RS model as a special case, with a scalar field 
in the bulk. In contrast to previous work we worked in the frame of the induced metric 
of one of the branes in order to study collisions, which are always singular in the Einstein 
Frame. After discussing the significance of 'supersymmetric' pairings of brane interaction 
and bulk scalar field potentials from a fully five-dimensional viewpoint, we derived an 
effective theory to discuss the low-energy motions of the system about a BPS ground state 
configuration with small SUSY-breaking potentials. It was found that, for this specific 
class of braneworlds, the presence of a bulk scalar field has little qualitative effect on the 
dynamics of the radion itself, i.e. the motion of the branes through the bulk. However, even 
in the absence of any perturbations around an FRW induced geometry, its contribution 
to the brane Friedmann equation causes the Hubble parameter to diverge at the collision 
even in the Brane Frame. By monitoring the size of certain functions of the field variables 
which are required to be small in the derivation of the effective action, we can identify the 
regions in which the MSA itself is valid. Whilst collisions are well-described by it in the 
RS case, the presence of the bulk scalar field can causes it to break down just before the 
collision. Although the RS model can be used to construct toy models whose regularity at 
the collision gives a link between pre- and post-Big Bang phases, where the Big Bang is 
identified with the moment of collision, this work implies that any attempt to generalise 
the model to include bulk scalar fields encounters singularities. Recent developments on 
modelling brane collision in the context of M-Theory [15] suggest that, with more care, 
these divergences might turn out to be removable; purely within the context of general 
relativity, however, they appear to be unavoidable. 
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7. Appendix 

Firstly, the Ricci scalar decomposes as 

R(g) = -8^ -12^ + 4 
a cr or 

where R is the Ricci scalar of the metric g^. Writing S = J d 4 x^—g C, we find 

'z 2 (x) 



R f Z2[ - x > 
£ g D — / dz a (z) 

K 5 Jzi(x) 



5 Jzi(x) 

l+l/3a 2 f l+l/3a 2 

-R (7.1) 



2k K j (1 + 3a 2 ) 

where £j = £(z = Zi(x)). The scalar field Lagrangian contains no x— derivatives (as shown 
in [16] there is no additional perturbation around the BPS field configuration that needs 
to be taken into account) and hence makes no overall contribution. 

The moduli kinetic terms come from the boundary action. For the z-brane the induced 
line element and outward-pointing normal are given by 

dsj = {afg^y + d^ZiduZi) dx^dx v (7.2) 
= hf v dx^dx u 

fdx - = (-I) 1 ; = {dz - d^dx") 

l + Vzf/a 2 



giYing V / ^)-a| v ^fl + -Lv^ (7.3) 



2«, 

where V is the covariant derivative of g^. Here, and from now on, we are assuming that 
Vzf <C of, i.e. that the brane fluctuation lengthscale is much larger than that of the bulk 
curvature. Then, from ( |3.12| ), 

L B D -of v - ok w + ^ (^alVzl - -?-a?Vz? ] (7.4) 

«5 \52 6. / 



where we have from (3/7) and (|3.9| ) that 
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and we have assumed that the SUSY-breaking potentials v and w are also small. Finally, we 
need to compute the Gibbons-Hawking boundary term involving the extrinsic curvatures 
of the two branes. Taking outward-pointing normals again, 



--d a 



-9 »fo 



Z=Zi(x) 



a 4 



i-iy 



lVz 



2 a 2 



-9<H 



St 



+ 



(-1)< 



fc- 



6 



/ / 2 = 2; (z) 



Z=Zi(:r) 



.)} 



Z = Zi(x) 



From (7.3) and (3.12) it then follows that 



£gh D — 



6k 



21 (a;) 



fl 2) 



A; 

af — Vz 2 - 2aiV^ziV^ai - (1^2) 



«2- " 



6 



(7.5) 



where integration by parts has been used in the second line, and using fl3.9| ). It is straight- 
forward to show that the omitted terms above some to zero as expected, leaving an effective 
action 

1 



2/l Hvt^ 



(fxJ-g WR + 



1 

6c? 



6 



l/3a 2 -l 



(7.6) 



f2 



1 



1 + 3a 2 



f l+l/3a 2 f l+l/3a 2 
SI ?2 



2/3« 2 tf 2/3a 2 



W 



As expected, the potential for the moduli fields vanishes for v = w = 0; when supersym- 
metry is unbroken (i.e. ( 2.1C| ) holds), the moduli fields are massless. Next we write the 
action in terms of the induced metric on the positive-tension brane, hfy 
order it is sufficient to use the approximation 



hfj, u . To this 
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Using 



R = a\ (R(h) - eaiDa^ 1 + ...) 



(7.7) 
(7.8) 
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we finally obtain, after integration by parts, the action as given by ( |3.14 ). 
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Figure 5: Divergence of r\ as ip — ► 0: For four different sets of parameters and initial conditions, 
77 can be seen to diverge at the collision with opposite sign to a. This divergence causes a corre- 
sponding divergence in the parameter e±, signalling a breakdown of the effective theory description. 
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Figure 6: Late time cosmological evolution: For suitable initial conditions ip is rapidly driven 
to unity and the branes move far apart. The scalar field grows logarithmically, the Hubble parameter 
is slowly driven to zero, and the Universe expands rapidly. Here a = —0.2 
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Figure 7: Late-time i](t): The evolution of —arj/2 as if> — > 1 is plotted, with a — 
0.4, —0.2, 0.1, —0.05 from top to bottom. The dashed fiducial line with unit gradient demonstrates 
the late-time behaviour 77 ~ — §logi. As a decreases it takes longer for the system to settle into 
this behaviour. 
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